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ABSTRACT 

We develop a theory of Orbits for the Inverse-square central force law which 
differs considerably from the usual deductive approach. In particular, we mate 
no explicit use of calculus. 6y beginning with qualitative aspects of solutions, 
we are led to a number of geometrically realizable pnystcal invariants of the 
orbits. Consequently most of our theorems rely only on simple geometrical 
relationships. Despite its simplicity* our planetary geometry is powerful 
enough to treat a wide range of perturbations with relative ease* Furthermore, 
without introducing any more machinery, we obtain full quantitative results. 
The paper concludes with suggestions for further research Into the geometry 
of planetary orbits. 



*This is a preprint of a paper to appear in The American Journal of Physics , 
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K SomcOualLtBtreBUcsults 



]. 1-nirodutfion 



From junior high school on, students of science are taught that K*pler h i Laws 
describe the motion of planets around the sun. They are given no hint of how they 
themselves art understand the "why" of these laws. By high school the students have been 
taught Newton's discovery, that the inverse-square force law accounts for those beautiful 
ellipses, but the connection is not yet for their eyes. After * year or so of college it's finally 
time to plow through the thoroughly standardiied and unmotivated proofs., using intricate 
manipulations with differential equations. 

In this paper we outline an approach to orbital mechanics which is, accessible to 
beginning physics students and presuppejrj no knowledge of calculus. We give an 
elementary (yet mathematically correct) treatment of Kepler's Laws and also investigate a 
simple first-order perturbation theory for orbits in an tn verse-square field. Our theorems 
and proofs arise naturally from trying to understand orbitj in terms of their physical 
invariants. We therefore feel that our treatment provides a better view of "what doing 
physics is really [ike" than does the standard route via algebraic manipulations. 

The key to the method lies in comdering the velocity Ipacc picture Tor a planet's 
motion about the sun. The concept of a velocity space is not normally encountered by the 
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student until he or sin* is presented with phase space in the context of the formalities of 
HamiStoniail mechanics. We think therr ii little niton for this delay because, though it is 
usually consjdered an "advanced concept", (he v*lKlty spies picture of a pamde's motion 
lies at the core of Newtonian Tmech an Lei. indeed, the qualitative content of Newton's F-ma 
{for otir purposes, miv-FAt) is simply lh;ir fihytital interactions between objects take place 
by a modification of velocity, ra:her than by a change in position, Appreciation of thli 
fact can greatly assist the development of physical intuition and understanding, and velocity 
Space Is a natural tool for exploring Newton's primarily conceptual breakthrough. In Part 
III, when we discuw perturbations, we shall jee what rich dividends can be yielded by 
looking at physical phenomena in the right conceptual frame-m this case, velocity jpa«. 



In the following presentation, we riH>e tried 10 waLk rather a narrow path between 
two extremes, On one hand, a description of our methods and results, would like no more 
than a few pages if we used the full precision or mathematical apparatus ^including 
calculus.) available to science Students after 4 few years of university education- On the 
Other hand, we could have spent considerably more space developing a complete and ieif j 
con taitied course for very early physio students. Since we I eel the material can be useful at 
both levels of physics education, we have attempted a compromise. We apologia both to 
those Who find our presentation extended and perhaps verbose, and to those who might 
find It sVecchy and incomplete. 



We gratefully acknowledge the inspiration arid encouragement of Seymour 
PaperL He introduced us to this way of thinking about orbits, and pointed out the basic 
results described in Section; 2 and 1 These sections closely follow pam cF his piper 
(reference !)> which preiwits * broader view of the conception of education in science and 
mathematics from which this work. grew. We would like to thank the editor and referee- 
From Tftt AmeHvi-Q- /puraail' af PkysiLi for many encouraging and helpful comments and 
also to thank Euiin jabari oF the M.LT\ Artificial Intelligence Laboratory for preparing the 
Illustrations for this paper. 

2. Thi Orbit U Cttoied 

Standard approaches begin with' the arduous task, of proving that planetary 

■ 

orbjts are precise ellipses. We begin by proving a more qualitative proposition* that orbiis 
are closed. In doing so we dispense with a great deal of analytic: clutter, and the important 

special nature of inverse-square orbits which nukes (hem closed comes into central rows, 

■ 

We will prove that no orbit like that in Figure I is possible. 

If a plantt crosses a Aa{f-Une ft&n tnt inn twite, tfun a ewjfj it ef 
iftt same p otn: cack timt^not further &xt &t doser in. 



We assume twp pieces of knowledge 




Figure 1: An Impossible orbit. 




Figure 2: Opposite pieces of orbit, at radii ^ and r £i sweep out areas 



A-, and A.£ in times At-, and At^, 



L The force on the planet,, when U is distance r from the surr, Is K'r^ towards the 



sun. 



2- Angular momentum j; conserved. We use this in the form of Kepler's Law 
that the radius From the sun to a planet sweep* out equal areas in equal times, This can be 

easiij 1 derived and we remind readers of iti simple geometric pronF in the appendix. 

Now consider diametrically opposite pieces of the orbi* which subtend the same 
(small) angle measured from the sun, as in Figure 2. Kepler's law states thai 



A 2 At 2 

What else do we know about the area or the time? Geometry tells us that 
A. 



i_ 
A 2 ' r = 

2 



Those r* r S arc ton suggestive for US not 10 make a connection with 
F = KJr ? , In face 



" =: — - — hence — — ■ — - — 



and we conclude FjAtj ■ F^At^. Since 7j and*F 2 pu[l in opposite directions, 
F, At, - -F 2 At 2 



W* can identiFy these terms: According to Newton's Second Law, on each piece f the 
orbit, F4t ii precisely At the change in vctod?^ which we call the "klc^ associated to 
that piece. Thus the tat equation sayj »hit the change in velocity over one pica? of orbit 
•KKlf, cancels the change on the opposite plea. Starting en the half -line which the orbit 
crosses twice, divide the orbit all the way anund into similar p* irS oF opposite pit™. The 
total change in the plants velocity tauten succeisive crossings of the half -line is the sum 
of the chants in velocity over each small piece; adding thwe up in opposing pairs, we see 
thai the total change IS icro . Whenever the plane! crosses a given half-iine h It has the same 
velocity. 

Now Kepler's dictum of equal area in equal time allows us to conclude that at rwu 
crouds D f the halHine, not only velocity distance from the sun is ihesarrte. Figure 3 
shows the areas swept out by the planet in snms short time it after successive crossings vf 
the half -line through A, B, and O. The Telccitiej at A and E are equal Therefnre the 
pieces of orbit AC and BD Hre both equal totAt, but the area oF AOC must equal the 
area of BOD. Then A equals R and the Drift doses. 



3. A Theorem ft! VelaUty SpOtt 

The preceding proof r«red m the fact that the kicks over opposite pieces, of the 

orbit have equat nmgrii'.udes; 

Fjrfl, = F 2 m, 




Figure 3: AC-flD, area Moarea BQD, therefore QA=0H and A-B, 



^ Atf 
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Figure 4: Part of an orbft divided into equal-angle slf 



ces 



But, to derive this equation it hn ' t hECS ^ ry thM fhe p ^ ^ ^ ^ ^^ ^ ^ 

only that 

2 



A r 



A, " r* ' 



and thil is [rue for.any two pieces of the whit over which the radiat angle changes by th* 
same mull amount. So, if W divide the whok orbit Lnto small p.eres subtending th, « me 
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angle A? p the "kick" v«loi Tor the various pi«« a f| f, ave the ume frngth (Figure 4). Not 
only a* alt the leitflh^ equal, bu; th, ration between suceisiv* kick v«tors U COmtatlS ind 

equals b$ . "' 

Thus w* have a very simple algorithm For generating th, changing velocity as 
the pkn* movts along iu orbit. Starting at a given velocity vector w add on kick vectors 
one After another. Each addition i, a step qF constant length and successive step* difter In 
direction by the constant turn, A* , k it Miv tD Ke lhat ^ alpiri[hrn . Grj FORWAED a 
short distance, TURN through a smiil angle. GO FORWARD the same short distance. 
TURN through the same small 4nf | e _ W J]| ^erare a circle. 3 We conclude that rh E kick 
vectors line up along a circle (Flgun? 5}. 

We car: intopret "adding on iU ««*ive kick vector" by introducing the notion of 
vtiacitj j^tf. The vetority of an object is usually described by a 'velocity teeter*, that is, * 
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Figure 5: Placing the kit* sectors end to end. 



direction and a speed {1^^. in comparing different velocities , c iB useP u | to pilt , hc , a „ 

of alt velocity v^torj do^n at some common point, 0, and to depict » velocity by the poi nl 

where the tip of the velocity vector lands. Wish this convention, we «n draw t Wo different 
picture* Ib describe the modern of an object 

JJ) the collection pf successive position* of [he object In >«r space; and 
K) the collection of successive positions of the tip of the object"* velocity vector. 
Thta is a path in "velocity ijJacc", a picture oF how velocity changes. 

The second picture is called the 'velocity jpaee pa th" or Velocity diagram 1 *. 
Figure 6 exemplifies thes? two kind* of diagrams. 

Velocity is the thing tint changes position, kick, are the things that change 
velocity. To get from an object's position at one instant, t( to its position at t-At, *r add c,r, 
tt>e vector vAt. To get from the object velocity at one instant to the velocity at a slightly 
later time, we add on a kick vector, Fit. Adding up successive vAt vectors give* the 
position space path: adding kick vector* give* the velocity space path. We can now restate' 
our result as 



Cttdt Thmem 
For an ctfecf ir.wmg tn an trntrnt-Hputn fitid r tin velxtij ffwtt 
path lies on. 4 ctrclt. 




VELOCITY PATH 



Figure 6: One object's position spare and velocity 



s puce paths, 



To avoid confusion we point out that the center of this circle K ml necessarily at 
the origin in velocity space. 



4- Thf Velocity Space Path 

Our Circle Theorem tells us that the velocity ipare path lies on a circle. But is. It 
a complete circle or just oak of it? We can answer that question atid a bit more. 

We showed that an orbit which does manage to get a]] the way around the tun 
crosses every half-line from the sun exactly once. Surh an orbit Is a simple closed curve. 
In a complete revolution, therefore, the direction of the planed velocity vector must change 
through 360*. (See reference 4.) That mean* that m vslataj ip&ta oto the path meets every 
hatHine from the velocity space origin. It follow that closed orbits in position space 
an-Bipond to complete circles in velocity space, and we have learned, besides, that the 
origin of velocity space J j inside the circle. 

■ 

We now have a good qualitative picture of L>ie velocity space, path for a closed 
orbit. [Open orbits are di mussed m section II.} In Part n, we will extract information about 

the position Space prfcut from our velocity space diagram. 
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Figure 7: D = radius x A & 



It Jnvartanfe ol ihe Orbit 

5. Angular Momentum 

Ai we have «en, the velotitjf spate path of a planet in an inverse-square field li E s 
tin S circle. One obvious invariant of a rirck ii ill radius. How can we interpret this 
Invariant physically? 

W> got the circte in wetlon J a& the rwult or rhe algorithm "forward distance D, 
turn angle AG, repeat." As one ran see (Figure 7}. thii generate a circle of radius BfM. In 
our case Afl was an arbitrary Small angle and D wai the magnitude of the kick FAe over 
the corresponding small piece of orbit Letting u denote the radius of the velocity circle, we 
have 



- Fit 
LI- 



AS 



It ii not immediately obvious, thai tbii ii a Constant. However, we tan Simplify the 
expression using the fact frnm geometry that th* area swept out over a small p]cee of orbit 
is 

Then u»E£i - TrM. We can eliminate the apparent dependence on the non-constant term 
r by HJlng F-Kjr^ we obtain 

K 



% 




The term 2AfAt, which telh ho* fait the planet if sweeping out area, is precisely the 

constant called anguiaf WvomejiSum, L. (See Append^,} Therefore 

The radius of the velocity circle equals (be force constant 
K divided by the angular momentum L: u^KfL 
For a fixed gravitational held, the radius oF the velocity circle lelli us the planet's angular 
momentum: a larger radius gives a smaller angular momentum. . . 

6. One niatlsn 

The velocity circle has another ^variant so obvious it Ji easy to DYflrfook-the 
position of its center in velocity space. As we remarked above, even though the origin j.j 
tnitde the velocity circle J t need not be at (As e^tar of the circle. Let t be the vector 
running from the origin in velocity space ^o the center of the- velocity circle, and let u be 
a radial vector oF the circle (Figure B).- Jn terms of tand u, we can think of the planet"; 
path in velocity space a? follows- at each moment the velocity v ii the sum of a comtent 
meter j and a vector u of anttant imgtA (*qu»l to u-K'L), v-l+u. The veic*iE ? space 
path ii generated as the radius "5. sweepi around the tip of the invariant vector^ 

There is a quite remarkable relaslon between the motion of f h the portion space 
radius vector (tall at the sun, head at the planet], and the motion of this "vekrclcj space 
radius." u. 

CarrtietioA ef Anglu in Potman anrf Velocity Spat< 
M each mvujtt thi plantfi radl&l wetor r u jvr frenrftcular to the 
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Figure 2 : The invariant "£. 



tadiui u oftht vtletlty fade. 
To see Chi* we examine how the Jack* fit into bath diagrams.- In position space each kick is 
para del to the radial vector r. in velocity space the same kick Js tangent to the velocity 
circle and hence perpendicular to the velocity circle radius u. 

Hence u is perpendicular to r. 

Correlation of Angles Itapowerful principle, It tells, use 

(1) Each, point on trie velocity space path corresponds to a unique point in the 
planet's, orbit. (The planet cannot attain the a me velocity at two different points in its 
orbit.) 

(2) The planet's position vector sweeps around the sun at the same rate and with 
the same direction (clockwise H counter clockwise} li its Velocity sweeps around the Circle 
in velocity spate. The two are always 90 degress out of phase. (Figure 9) 

Now we. can.gjve more meaning to the z sector. The planet's speed, v, is 
greatest when the u and i vectors are lined up, least when they are opposed; v i -u-z, 
V min' Ut1, Therefore the t vector points in the direction of maximum Speed and 
□ppo-site to the direction of minimum speed, 

The points of .greatest and least speed occur where the velocity vector is parallel 
to u, and SO perpendicular tp the position space radius"?". It is not hard to show then 
that the point where Speed attains its maximum (respectively, its minimum) corresponds to 
the minimum (respectively, maximum) distance from the sun. The reader can fill in the; 
details of the proof sketched below: 
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VELOCITY 



Figure 9: Snapshots of position and velocity over an orbit ujj. 



Step 1: At a maximum fir minimum distance the velocity y mutt bt 

.1 
perpendicular to the rarijws r. 

Step 2. The velocity diagram shows (hat there are precisely two points where this 
ran occur. 

Step 3. Conservation of Angular Momentum implies that r maJ( corresponds to 

Y min * nd r rnin to W' 

Now we- have same more qualitative information about the shape of the orbit: 
there is precisely one point of maximum distance from the iun, and one of minimum 
distance, They occur on ori^os-itc sides of the sun since the corresponding U Vectors point 
in opposite directions. The~i Vector determines the orientation of the orbit. It points in 
the direction of maximum speed, (Figure 10) 



7, Shape srtct Symraifry 

What more does the length pf i tell us about the orbit? Consider what would 
happen if i vanished. There would then be no direction picked out for maximum speed. 
or distance from the sun. The planet wpuld hive to travel around the sun In a circle at . 
uniform speed. (Another way to see this: % would be ectua! to lT and therefore always; 
perpendicular tol 1 and of constant kngth-trte characteristic of uniform circular motion.) 

This suggests that i indicates how the orbit deviates from a drcle. We can 
.make thii precise. One obvious measure of the non-circularity oF the orbit is the diFFerentt: 
in the extreme distances fro^j the sun, 




max 





Figure TO; t determines the orientation of the Orbit. 



r maK~ r min" 
If we want an invariant iba'. depends only on the shape and nw the sue of the orbit it is 
better En see how much the ratio 

r maj/ r mm 
diFfers from t. The lerifS.ru oF i, relates ihr maximum and minimum Speeds; 

T min- uz W 1 "* 

To relate speed to distance from the sun we use angular momentum. If the angular 

momentum is L then 

mm 'mat man * r min 
because at these places in the orbit v is perpendicular io r (Section 6) r Therefore: 

u+z 
u - z 



"max 


Vmiix 


1 min 


Vm in 


Since u-K/JL we can 


also Write 


^max 


K+Lz 



But K depends, cwily an the nature of the gravitational f Jeld ao we sec that nur "shape 

invariant" is determined by Li, The larger Li, the more the orbit deviates irum i circle. 

In section II vte will derive the analytic result that the oibit is an ellipse and Lz determines 

its eccentricity. 

In fact., knowledge oF K and the velocity diagram essentially determines the 
. planet's mrr.ion. The maximum verity can be read off the diagram immediate^, as can 
, U-Kr'Li'So we fcnpw L. The shortest radius in position spate has length Thi^L/Y^ and 



Is perpendicular to t mas (Section 6). Having this one vector,"^, we can generate the 
orbit starting at the position determined bjr'f |nLn with the foltaiwjng algorithm. 
E. Travel a shot! distance vAl in the direction of V, 

2. Measure the change in angle, M, m position space 

3. Firtd the velocitj it this new angle {by rotating 'u I h rough A* and consulting 

the velocity diagram), 

1 Return to step |. This generate! the entire- posrtion space path, 

■ ■ 

Notice (hat the velocity diagram is symmetric about the line determined by z. 
The above algorithm translates this tact into a symmetry of the position space orbit. 
Starting at the nearer p Jin t tn the *un. conduct the orbit in the forward direction for a 
while, along vj for M [h then along \ for A^ h and so on, Now go back to the 
starting point and run the algorithm bacfc wards with the nine seqtieocs of M'j, Sinte the 
Velocity diagram ll symmetric we generate the saraesrcift" segments of prblt, except that they 
hive been flipped about the line perpendicular tot Therefore the entire ert.it is 
symmetric abouic this line, . 



[ram. 



We have so far obtained the following information from the Velocity dJagi 

1. The radius of the velocity circle determines the orbit'* angular momenium: 
u-K/L, 

2. The center of the velocity circle determines the orientation of the orbit (f 



points in the direction of maximum speed) and its "jTiape* <Li determine* r^^r s ). 

3. From Che velocity diagram, ws can algarilhrnjcalty determine Lhc whols orbit. 



IN. Perturbation 5 

9. Tht PtTtuxHUdn Formula; Re4(^ TftmiC 

It h in the study of perturbations, or how orbits change under small kicks other 
than thtwe given by the sun, chat our use gf veiodiy diagram* really pays off. It pays off 
for a very good reason. which we mentioned in the introduction as a qualitative- form of 
Newton's Second Law of Motion: 

Force acts on trie path* of parLLcJes by changing vtfauy and not 
poititm. 

If we fail to take account of this fact we may be faced with situations that appear 
BounDwnntultlvfc For example, suppose a spaceship in a circular orbit around a planet 
applies a small outward thrust (Figure llajf ' 

How will the orbit change? "Intuition" may suggest that the orbit will elongate in the 
direction of the thrust, something like Figure lib. In fact, tne orbit will elongate, bur in a 
direction perpendicular to the kick a? in Figure lie. 

To understand this we consider how the kick changes the velocity diagram. The 
space-ship started in a circular orbit whose velocity diagram is centered on the origin. 
Since force affects velocity and not position, it is reasonable (and we shall show below) that 
the affect or the kick in velocity space u&Uy u just to move tne velocity circle in the 
clirecclofi of the kick (Figure 12). The corresponding change in the pwltion space orbit is 
the "counter-intuitive" effect described above. 

Our strategy for studying perturbations will be to see how kicks change the 
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Figure 11: Will the outward kick dn orMt (a) produce (b) or (c)7 



AAv 





POSITION BEFORE 



POSITION AFTEfl 




^' ■ ji^l 


-1^ 




/ 1 * 


■^ 

-v 


\ 


/ s 


M 


1 / 


\\ 


// 


\\ 


1' \ 






\ 






l\ 


l\ 






_— -r K 


^U ; 


7* "~ 




■" 


. 







VELOCITY BEFORE 



VELOCITY AFTER 



Figure 12: The perturbation induced by an outward fcfck. 



TdflCity diagram More precisely, we know that the shape and orientation of the orbtt is 
determined by Li the i vector time* [lie angular momentum, so we want to Find the change 
ill Ll r A(Lz), produced by an arbitrary kick. 

The bailc velocity space Equation, v - 1 * u, gives Li - L v - Lu. Since u naj 
l*u£th EUL we have Lu- Ks, where tit a vector or unit length whose direction Is 
determined solely by the abject 1 * position. (It is perpendicular to the radial vector-*) To 
compute the effect of a kick At on Lf - Lt . U notice that s.n« kic^s do not affect 
position, 1 is unchanged. K is also unchanged. Therefore the changed lA H the same as 
the changs in Lv> and the first-order approximation to the change (ft a product of 
changing quantities gives 

■ 

Penttrbatlsn Formula: Ai^-'ti&l;* LA^v 

We can use the Perturbation Formula to tidy upaur discussion of the "radial 
thrust problem* (Figure 1ft Since Lv Is in the radtal direction, the angular momentum 
does noi change (AL - 0), so the formula implies A(Ll> - Li v. This means that the 
velocity diagram change! from at Vector of zero Co at vector tn the direction of A"? 
{Figure 12). 

Intuitively, th? velocity circle Ji "pushed" in -he direction of the kick. Note that an inward 
kick at the bottom of the position spite orbit would have the Same effect. 
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Figure T3: The perturbation induced hy 4 tangential kick. 



VXTangtmal Tfrust; Solar Wind; t^^Lamt 

In this fatten we apply the Perturbation Formula to some other orbit problems. 

Tangcmiat r^uj/; Suppose again that a rocket starts in a. circular orbit, but this 
tin* provides t tangential kick (Figure 13>. To determine A(Li) - ^Al ♦ LAv" we note flat 
vAL js an impulse in [he direction of Atf since AL is positive and v is parallel to A y r 
Then the newljr crfttH3 £ VcctQr mUs[ he jn thp d ^ ctjOR rf lhp ]mpube ^^ ^ ^ 

elongation in position ^pate j s again perpendicular to the kleL 

r*< 5oTflr ffJnd? Assume the. recta is affoote4 nor only by the planet's gravity 
but ato by a small constant force (Luehrwm's 'solar *ind"). If the perturbing Force is 
small compared to gravity, each revolution or the rocket will be nearly an tltipw. W* «n 
therefore think of Lh* orbit u an etlipje which varies through time. To compute how the 
ellipse changes we view the wind as providing impulses all aloog the orbit {Figure I4> and 
sum A(Li) -. vAL + LAvover one re volution. 

The LA* contribution is a net change in the direction of tt. To compute vAl 
we notice that AL is positive on the bottom halF of the orbit and negative on the top half 
u shown in Figure 15b. We can sum the tAlA by exploiting: the symmetry of the orbit, 
The vertical component* of the vAl/s on the feft cancel the vertical components of the 
vAL'i on the right, leaving oniy a horjwntif component in the direction of At {Figure 
l£c>. This adds with LAv to produce a t vector in the direction of the wind. fntuitively, 
the velocity circle gets *bl&W in the direction of the u,inc3. The orbit elongate, 
perpendicular to the wind as m Figure |$, 
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figure 14: The solar wind. 
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F^ure 1&j The Vectors *v (a) and t^L (h> for the solar wind. In fc) 
we see that the vertical components of 'vAL on the left 
cancel the vertf cal corrpOTients on the right. 
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Figure 16: Position space change in the orbit under the solar wind. 



Since the wind does not afFect the symmetry of the orbit about trie vertical a^Jl we can 

I ■*_ — % 

apply the same analysis as above to show thai the Li vector continues increasing in the 

direction of the wind. From this we conclude that the orbit becomes more- and more 
eccentric while the direction of *i remains COhslant, and the orbit continues to elongate 
perpendicular to the wind, 

The orbit becomes closer and closer to a straight line, and it eventually reaches a 
point where the "snuff wind can have 3arge qualitative effects over a timescale of less than 
one revolution (the orbit in fact reverses direction), and out 1 method of averting over an 
entire revolution becomes inappropriate. 

Tie r -t Paw Fit!d: If * is a small constant (we will take it to be positive), we 
hi treat the central funce field of magnitude r"^ +E ' at a perturbation of the f^ field.. The 
perturbing force is some force (positive or -negative)' in the radial direction. To understand 
how this perturbation affrcts the orbit,, we make the important observation that the shape 
of a r" or r~W**f orbit does not depend on th« scale which we use to measure radius,, 
Therefore we can determine shape by using any scale which males it convenient to 
cotnpuite the effect of the perturbation. For the orbit shown in Figure 17, we scale to make 
the distance OP equal to one, Since 1/r 2 < l f r^ +(3 for r < I and l/r 2 > ti/ 2 *^ for r > I the 
perturbing force is as shown. 

For thii perturbation the kicks are radial, so L is constant. This means A(ut) - 
LA£ but from the perturbation formula, AL - implies A(Ll) - LA£ Hence Lt ~ AvT 
Now we can sum At -Av^cver an entire orbit. The leFt-righi symmetry of the crbit and 
perturbing force means that the sum of horizontal components of the kicks must cancel 





Figure ]?; Orbit with perturbing force indicated for r 
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Figure 18: &~z perpendicular to z. 




Figure 19: Precession , of r"' ?+t ' orbit. 



I it doW nward, that is, perpendicular to the original ;.. Subsequent ATs wilt be 
perpendicular to the current Tand this results primarily in a rotation off. not a charge Jn . 
length (Figure IS). The 'major axis" of the position ipace orbit, by the consequences we 
derived from Correlation of Angle, must follow this counter-clockwise rotation. Though 
the orbit retains it; shape Li r it processes (Figure 19), 

Warning rlt should be remarked that in the preceding two examples we tooled: at 
the F vectors as representing iv for the perturbation formula. Of course we should hive 
used Fit r but, because of the iymmetry involved, the At factor can be ignored in those two 
cases. In more complicated situations, though, this dops become an issue. For example, we 
incite the reader to use the technique of this section to treat the perturbation Induced bf" 
an ablate iunr 

TV. Analytic HesuLrs 



II. Tht Orhti ts a Cmtt Statpn 

An objection that is sure tooctur lo some of our readers: goes something like thU; 
"All these intuitive methods are fine, but if )'ou want useful quantitative in Format ton you 
have U? return Co the itandard differential equations youVe been trying to get along 

with nut" 



Of cou:« there are orbit problems our simple methods won + t handle. As for the 
itandard results, however, we arc able id derive the orbital equation directly from nur 

velocity diagram using no more than tngon-Ometrj 1 : 

The orbit is described In polar Coordinates by the equation 

L 
r ~ — -^— ^~ 

The proof IS a natural correlation of the basic quantLtkJ t r h ~v„ u d I and L tiling 

the definition of angular momentum At a point th ibe o*-bit when the"y vector and the~z 

vector differ in direction by an angle S we construct the angular momentum triangle (see 

appendix). 



The area of the triangle OAB in Figure 2D is t>j definition Lf2. If h is the height 
of ihe triangle t>Ffen 

1 . i 
r L= _ rn 

2 2 

Since u. and r are perpendicular, the height or the triangle U given by b - u - t wi #, 
Therefore 

L - rh - r {u - z cos ft 

Here fl represents the angle In vgf.K£ty J/#« between U and Lhe fixed vector t Correlation 
of Angles implies that & alio measures the angle in position space from r to a fiKed ¥*ctnr 




~z cos9 



Figure 20: The angular momentum triangle OAB. 
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Figure 21 : Sample arbfts in position and velocity space. 



perpendicular to t. Therefore r and art potar coordinates Irr position space. 

The abovE equation describes, a conic section. When the origin of velocity space 
lies within the circle, u>i and thecrhit ti in ellipse. When the origin j* outside the circle, u 
< z and the orbit is hypwbatk. When the circle passes through the origin, u - t and the 
orbit ii parabolic {Figure 21), 

Writing the equation in the form 






_ % 



l /k 



-(*) 



eos0 = 1 - 



( L %H* 



Lz 



we get Che "stand and form-" for a ennic section and see that / K is the eccentricity K and 

■> 
L /K 1j the radius of the orbit when the eccectricltjc Is. ienjL 



12. C<m titi/auoTi of Energy 

Energy ccm Serration d«i no( arise na:ni.i II y using this gEome^ric approach 
although we can obtain the result as a ilmple application, 



Apply the law of WJlntl tn the velocity diagram in Figure 22 to get 




Ffgume 22i v 2 ± Z Z + u z _ 2liz cosA 



3 _ 2 2 

v i +U — 2uz cos 6 



L 
Substituting i cosfl - u - y (from Section 11} we obtain 
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and hence (since u * 


K/L) 
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Since i and u are constant fo? the orbit, 2 is a constant, the total energy, E. It is 

interesting tq pole that when the planet trasses the wmi-mtlW aKts (F perpendicular to'v)"u 1 
X und "v form a right triangle with v 2 -u ? -^, hence the kinetic eneigy /2 is exactly the 
negatLvE of the tocal energy there. 

13. Kepler's Third Lam 



Vft can Use the relation of angular momentum to arEa swEpt out, 2A-Lt h to 
compute the period of the planet 4 * revolution. ]n one period the planet sweeps out the 

■ 

entire are* of its elliptical orbit The area of an ellipse of it ml -major aittj a and 
eccentricit j" e is given by A - r a *\ 1 ™B For I h e orbit we ha v e 
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Then [he period Ji determined ty 

LT= 2*V.Jl-e a ^ 2 * a 3 4 /JT 



L -J'.i 



or 

3 



T = 2^l!l^r - 21 % 



■a 



In terms of quantities appearing in [he velncify diagram we get 






li, t?p#rt Or*«j " 

■ 

r 

Rather than treat the hyperbolic case in detail, we leave the reader to verify the 

following; 

t For an open orbit, the arc of the velocity circle which it actually traversed Js the 
part shown below, bounded hy the tangents io the circle through the origin of velocity 

*p*«. Velocity Spate geonieEr ? ^vts the correlation of energ? with limiting velocity v 




Figure 23: The velocity diagram for an open *rbit< 




Figure H\ Deflection angle for an opert orbit, 



-jz 2 -u 2 -J 2 E~(Figure 23). 



2. The deflection angle (angle between the two asymptote* of the hyperbola} cm 
be easily found as a function of energy and angular momentum: ran — ^t^H (Figure 

3 K 



B, Sllggeltl&til for Furtfar RaartA 

We have by no means ejihauited the study of the geometry of orbit* in thli 
paper. The geometry of orbits, particularly the perturbation theory, Is a rich source of 
problems, even of mini-research projects of the type described! by LuehrmmrA Below we 
make 1 some suggestions ffir problems and Study tapirs. 

An tniiruttlvt ptiTadux.Th? Galilean transformation requires the relation between velocities 
measured m two different frames of rcfiente moving with relative velocity % to be"v' -~v 
* v , This it h a simple matter to mqve into a frame with relative velocity -t and transform 
the z vector for an elliptical orbit to ft Why then does one not observe a circular orbit in 
the new frame? In particular, what fails in the a1 E orithm pf Section 8 which don generate 
ft circle in position space givrn a circle centered about the origin In velocity space? 
Ah Aid to AsirogttfioT!? Suppose we had to pilot a spaceship in a gravitational field (such as 
simulated in computer "space war" games). Would a velocity diagram be a useful addition 
to oor instrument panel? For exampFe, to change from an elliptical orbit tp a circular orbit 



we fioed only consult the velocity diagram and apply a fare* to cancel the I vector. On the 
other hand, lots of information Ji lacking if i¥c use tatty Ihe Velocity picture- Intercepting 
another spaceship is a tricky problem invoking timing. (Although merely matching its 
orbit ii easy.) What other instruments should supplerMnt or possibly replace a velocity 
diagram'? 

GrttHtrtty of ikt kn-sRuitic sicillatw; \ fundamental geometric property of solutions, to the 
\fi dirrerential equation is that they have a vector con^artt of motion (^or (he maximum 
vebocity, or the Runge'Leni vector are ail possible choices for this, constant). The two 
dimensional harmonic oscillator with equal mode frequencies has a similar structure. 
Solutions have an obvious a a; is which may also be assigned a magnitude in any number of 
ways. Can one develop a useful velocity space geometry for that system? Can one treat 
simple perturbations, as with orbits? <*"* 

Mart $ trier Wind.- A further discussion of the solar wind phenomenon could mate use of ■ 
the fact that the force field is conservative and therefore the energy of the orbit is constant,* 
This implies that the angular momentum decreases as the orbit becomes more eccentric 
Thus the velocity circle is not only "blown by the wind" but also the radius becomes larger 
and finally infinite when the orbit degenerates to a line- Show that the -changing 1 velocity 
circle always passes through two fixed points in the plane. (Figure £51 

Our method of averaging over entire orbits is only a first-order perturbation 
theory whereas the formula 

d(Lz) d^ . dL 
— ■ - ' L — ; -+ V 
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Fijure 25s Changing velocity circle utirfer solar wf 



nd. 



fe is any parameter) is exact Orre would flke to have a treatment tf the w]ar witl d which 
Works near the turnaround prints, ifu approaching I, 

FinaElj, is there a complete perturbation theory based Pn the geometry of orbits? 
In particular Tiomk- can one treat perturbations out of the plane of the orbit? 



ffppgridix: Angular Mentntum and Ktpini Smmd Uu 

Throughout this paper we have b«n assuming Kepler's Second Law. There is a 
simple geometric proof of this which we place here in in appendl!t because ft did not 

originate with us. It «n be fqurtd in Newton's Frmtipla! 

The angular momentum, which we denote by L, is defined to be twice the area of 
the triangle determined by the velocity vector and the radius vector from the sun to Eh* 

planet. Ai shown below L is constant if the velocity doefl^t change - the triangles have 
equal areas since thev have equal bases (the length of v) and eoual heights (Figure Ml 
More remarkably, L T&mtns wmtant tf w changa tke velmty by applying <™ 7 kick m tht 
radial dimum (leuarfa & amy frm ihcwi). The a ff« t f a kicfc Av on the anguh 
momentum triangle is illustrated below (Figure A2>. The kick changes v to^v'. but Ermigl* 
OP Q, and OPQ! have the same oase, OP, and the same height (h in the diagram) since OP 
And ^C L J .1-, p ara i: Pl ThcrefnreOP^ancI OFq, f h,vv, the sam* i.rea, and ^rrgular 
momentum is unchanged, 

A planet moving about the sun, subject io no force but the sun's gravitation, has 
every kick applied in the radial direction. Nnne or these change the angular momentum, 
which Is therefore an invariant of the planets orbit. To find a geometric interpretation of 
this face, we examine the orbit- at time intervals it small enough that the velocity does no* 
change much over each Interval In each interva* the radius vector sweeps out a small 
triangle. The area of one of these small triangles is 



in 



lcs. 




L=vh 
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? 



figure Al : Triangles o£ equal area 




Figure A2 ■ Area, OPQ » Area OP.Q 1 
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(h as in Figure AIJ. The torsi area, s^tpt out over sornr long- rime 

■ 

U the sum of the ixeai of the &rna11 triangle* 



T" i i i 



1 1 

A= -Lit. *-— L A t 
" — LT 



Thi* gives Kepler's Second Law 

For a body moving in a radial fanoe ricld.. the- radial vector sweeps 

CUE equal area* In equal time^ 
It is unfortunate thai this prcwF is not more often presented in physio courses (although 
Feynman disojiiR ft and there is a movie 7 demonstrating this argument). 
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